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We analyze the Sudbury Neutrino Observatory (SNO) and Super- 
Kamiokande (SK) data on charged current (CC), neutral current (NC) and 
neutrino electron elastic scattering (ES) reactions to constrain the leading 
weak axial two-body current parameterized by Li^a- This two-body current 
is the dominant uncertainty of every low energy weak interaction deuteron 
breakup process, including SNO's CC and NC reactions. Our method shows 
that the theoretical inputs to SNO's determination of the CC and NC fluxes 
can be self-calibrated, be calibrated by SK, or be calibrated by reactor data. 
The only assumption made is that the total flux of active neutrinos has the 
standard spectral shape (but distortions in the electron neutrino spec- 
trum are allowed). We show that SNO's conclusion about the inconsistency 
of the no-flavor-conversion hypothesis does not contain significant theoretical 
uncertainty, and we determine the magnitude of the active solar neutrino flux. 
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I. INTRODUCTION 



Recent conclusive results from the Sudbury Neutrino Observatory (SNO) have estab- 
lished the existence of non-electron active neutrino components in the solar neutrino flux 
[1] and hence have given a strong evidence for neutrino oscillation. These results are based 
on the three reactions measured by SNO to detect the solar flux 

u^ + d p + p + e- (CC), 

u^ + d^p + n + u^ (NC), (1) 
i^x + e-^ Ux + e- (ES). 

The charged current reaction (CC) is sensitive exclusively to electron-type neutrinos, while 
the neutral current reaction (NC) is equally sensitive to all active neutrino flavors {x — 
e, /J,, t). The elastic scattering reaction (ES) is sensitive to all active flavors as well, but with 
reduced sensitivity to z/^ and Ur- Detection of these three reactions allows SNO to determine 
the electron and non-electron active neutrino components of the solar flux, and it is then 
obvious that the cross sections for these three reactions are important inputs for SNO. The 
cross sections for all three reactions are determined from theory, but the CC and NC cross 
sections involve nuclear-physics complexities not present in the ES interaction description. 
Thus the CC and NC cross sections have become the main source of theoretical uncertainties 
for SNO. 

The complexities in the CC and NC processes are due to two-body currents which are 
interactions involving two nucleons and external leptonic currents. In the potential model 
approach, the two-body currents are associated with the meson exchange currents and can 
be calculated in terms of unknown weak couplings. In effective field theory (EFT) , the two- 
body currents are parameterized. In both cases, experimental data from some other processes 
are required in order to calibrate the unknowns in the problem. In EFT, this calibration 
procedure can be described in an economic and systematic way. The reason is that, up 
to next-to-next-to-leading order (NNLO) in EFT, all low-energy weak interaction deuteron 
breakup processes depend on a common isovector axial two-body current, parameterized by 
Li^A [2] (sec more explanations in the next section). This implies that a measurement of any 
one of the breakup processes could be used to fix Li^a- A summary of the previous efforts 
in the determination of Li a can be found in Ref. [3]. 

In this paper, after briefly reviewing the EFT approach, we will present the constraint 
on Li^A using a combined analysis of the CC, NC and ES data from SNO and Super- 
Kamiokande (SK). Wc then compare this new result with other determinations of Li^a and 
comment on the interpretation of SNO's measurements with the assumption about the size 
of Li^A eliminated. 

II. EFFECTIVE FIELD THEORY 

For the deuteron breakup processes used to detect solar neutrinos, where the neutrino 
energies E,^ < 15 MeV, the typical momentum scales in the problem are much smaller 
than the pion mass m^(~ 140 MeV). In these systems pions do not need to be treated as 
dynamical particles since they only propagate over distances ~ 1/m^, much shorter than 
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the scale set by the typical momentum of the problem. Thus the pionless nuclear effective 
field theory, EFT(7f) [4-9], is applicable. 

In EFT (;("), the dynamical degrees of freedom are nucleons and non-hadronic external 
currents. Massive hadronic excitations such as pions and the delta resonance are not dy- 
namical. Their contributions are encoded in the contact interactions between nucleons. 
Nucleon-nucleon interactions are calculated perturbatively with the small expansion param- 
eter 

Q.<i^ (2) 

which is the ratio of the light to heavy scales. The light scales include the inverse S-wave 
nucleon-nucleon scattering length 12 MeV) in the ^5*0 channel, the dcutcron binding 

momentum 7(= 45.7 MeV) in the ^5*1 channel, and the typical nucleon momentum p in 
the center-of-mass frame. The heavy scale A is set by the pion mass m^. This formalism 
has been applied successfully to many processes involving the deuteron [9,10], including 
Compton scattering [11,12], np — > d'j for big-bang nucleosynthesis [13,14], i/d reactions for 
SNO physics [2], the solar pp fusion process [15,16], and parity violating observables [17] . 
In addition, this formalism has been applied successfully to three-nucleon systems [18]. It 
has revealed highly non-trivial renormalizations associated with three body forces in the Si/2 
channel (e.g., ^He and the triton). For other channels, precision calculations were carried 
out to higher orders [8]. 

For low energy deuteron breakup processes, it is well known that the dominant contri- 
butions to the hadronic matrix elements are the ^Sq transitions through the isovector 
axial couplings. The ^Si state (such as a deuteron) has spin S = 1 and isospin / = 0, while 
the ^Sq state has S — and I — 1. Amongst the spin-isospin operators 1, r", cr* and r^cr', 
only the isovector axial couphng r^cr* can connect ^Si to ^5*0 states. The ^Si — > ^5"! transi- 
tions are suppressed at low energies because i) the isovector operators do not contribute (the 
transition is isoscalar) and ii) the matrix elements of the one-body isoscalar operators vanish 
in the zero recoil limit {d and np states are orthogonal in this limit). This leads to large 
suppression of the isoscalar two-body contributions through the interference terms. Also, at 
low energies, the non-derivative operators are more important than the derivative operators. 
Thus the leading two-body current contributions for low energy weak interaction deuteron 
breakup processes only depend on a non- derivative, isovector axial two-body current, Li^a- 

In Ref. [2], EFT(;f) is applied to compute the cross-sections for four channels (CC, NC, 
Ue + d — > e'^+n+n and V^ + d V^ + n+p) to NNLO, up to 20 MeV (anti)neutrino energies. 
As already mentioned, these processes have been shown to depend on only one parameter, 
Li^A- This dependence is subject to an intrinsic uncertainty in our EFT calculation at NNLO 
of less than 3%. Through varying Li^Ai the potential model results of Rcfs. [20] and [21] are 
reproduced to high accuracy for all four channels. This confirms that the ~ 5% difference 
between Refs. [20] and [21] is due largely to different assumptions made about short distance 
physics. 

The same two-body current Li^a also contributes to the proton-proton fusion process 
p + p ^ (i + e^ + z/g. This is the primary reaction in the pp chain of nuclear reactions that 
power the sun, reactions which in turn generate the neutrino flux to be observed by SNO. 
The calculations in EFT(7f) were carried out initially to second order [15], and then to fifth 



3 



order [16]. Thus a calibration to SNO's CC and NC reactions can also be used to calibrate 
the proton-proton fusion process. 



III. FIXING Li^A FROM A COMBINED NC, CC AND ES ANALYSIS 

In this section we present the constraint on Li^a obtained from a combined analysis of the 
solar neutrino fluxes measured by CC, NC, and ES reactions. In SNO's analysis, a specific 
Li^A was chosen in CC and NC reactions. The extracted solar neutrino fluxes from CC and 
NC were then compared to each other and to ES to extract a consistent set of neutrino 
flavor-conversion probabilities and to map allowed regions in a 2-mass mixing description. 
Here we take Li^a as a free parameter and use the available experimental data from SNO and 
SK to fix not only the flavor-conversion probabilities but also Li a- The only assumption 
we will make is that the total flux for the active solar neutrinos has the standard shape. 

In the two-flavor oscillation analysis there are three parameters extracted, Am^gi ^12, 
and ^i,^, which are the differences between the squares of the neutrino masses, the mixing 
angle, and the total active neutrino flux. There are three separate experimental inputs, 
the ES, CC, and NC rates. It might at first be thought impossible to extract a fourth 
parameter, Li a, without additional inputs or assumptions, such as fixing the shape of the 
electron-neutrino spectrum. The shape is experimentally determined, but not yet with high 
accuracy. Our strategy is to note that, in active-only oscillations, there is no shape distortion 
in the total flux, and that the integrated spectral response in the CC reaction over a certain 
range of final electron energies is the same as that in the ES reaction over a different range 
of energies, independent of distortions of the neutrino spectrum [23] . 

The CC and NC are measured at SNO and the ES is both measured at SNO and SK. 
The measured event rates are the integrals of the effective cross sections weighted by the 
solar neutrino fluxes that reached the target. 



where Ri is the event rate, /j/. is the Ui flux and E is the neutrino energy. 5j is the effective 
cross section, defined in Appendix A, with i — e, ii,t for i'e,ij,,T + e ES interaction. These 
effective cross sections are the true cross sections convoluted with the detector resolution 
functions which describe how the energy is transferred to electrons and detected by their 
Cherenkov radiations. The effective cross sections depend on the electron detection threshold 
rpth_ For CC and NC reactions, they also depend on Li^a- 

The total fiux for the active solar neutrinos is assumed to have the standard shape. 



where 08^ is the normalized shape function (J^ '"'^ dE(j)^B (E) — 1) [22] and is the 
magnitude of the 1/^ fiux. This assumption is valid if there are no oscillations to sterile 




(4) 



(3) 



(5) 



{E) = $,,085 {E) , 



(6) 
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neutrinos, or, even if such mixing is present, the survival probabihty to active neutrinos is 
energy independent. Similarly, the flux is 

f^AE) ^ ^.JsB (E) P^^/^JE) , (7) 

where Pue/v,: (E) is the probability distribution of finding a out of a i^x- Obviously, 
Pu^/ux {E) is bounded between and 1. 

We follow Villante et al. in Ref. [23] to define the averaged effective cross sections Wi and 
the normalized response functions pi (E) for the spectrum 

dEui(j)&B {E) , 

p.{E)= . (8) 



/ 



dEai(l)SB (E) 



Then we rewrite eqs. (3-5) using eqs.(6-8) 

Rnc = '^ux'^NC , (9) 
Rcc = ^^x^cc J dEpccPu,/^^ , (10) 

Res = + J dE [aePe - ^i,,tPi,,t] Pu^/u^ ■ (11) 

An important observation made in Ref. [23] is that the terms with the normalized response 
function p (E) dependence can be related by choosing suitable detection thresholds. This 
allows us to reduce the numbers of unknowns such that eqs. (9-11) become solvable. These 
approximate relations and their corrections are systematically explored in Appendix A. As 
shown in Fig. 1(b), pcc is very insensitive to Li^a- Thus we can set 

d 

Pcc {E) |r*,''^=5MeV = . (12) 



Also, as shown in Fig. 1(c), to a very good approximation, 

Pe {E) |T«''=6.8MeV — Pfi,T (E) |rt^^=6.8MeV ■ (13) 

The corrections of the above relations change Li^a by up to a neghgible amount of 0.25 fm^. 
A more significant correction comes from 

P (E) = Pcc (E) |T*/j,=5MeV " Pe {E) |r|''=6.8MeV • (14) 

We will use a parameter e in eq.(19) to parametrize the correction. 

Because of eq.(12), the Li a dependence only shows up in ajsic and acc- The scaling can 
be written as 

'^NC — o^nc9nc{Li^a) , 

o'cc = ^c9cc{Li,a) , (15) 
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where the and o^<c are the values used by SNO [1] which are based on the calculation of 
Ref. [24] with the electromagnetic radiative corrections of Ref. [25] (see Ref. [26] for an earlier 
attempt) and a 5 MeV electron detection threshold. These cross sections are corresponding 
to the NNLO EFT results with i>i^^(7I = m^) = 4.0 fm^, where the renormahzation scale /J 
is set to the pion mass. In the following expressions, we will suppress the JI dependence of 
Li^A for simplicity. The scaling functions can be parametrized as 



9nc{Li,a) = 1 + aivc - 4.0 

9cc{Li,a) = l + acc - 4.0^ , 



(16) 



where ajvc = 0.013 and acc = 0.010. It is interesting to note that if a were averaged 
true cross sections instead of the averaged effective cross sections, then ajvc and acc would 
be almost identical [2]. The larger difference here is due to the difference in the detection 
methods. As shown in eqs.(A2) and (A5), the CC event detection requires the final state 
electron energy to be above a certain detection threshold, thus leptons transferring too 
much energy to the hadrons will not be detected. For NC detection, however, there is 
no such discrimination — all the neutrons generated in NC have the same probability to be 
detected in the thermalization and capture process. Thus in general, the scaling factor 
of NC is associated with that of more energetic scattering (with larger energy transfer to 
the hadrons) than that of CC. We find the qualitative difference between a^c and acc is 
consistent with the neutrino energy dependence of the scaling factors calculated in Ref. [2] . 
Substituting eqs. (12-15) into the eqs. (9-11), we have 

Rnc = ^i^x^nc9nc{Li,a) , (17) 
Rcc = ^.ycc9cc(L,,A) J dEpP^j,^ , (18) 



Res = ^v^(7e 



^"^+(1-^) (1 + ^) jdEpP.^i.^ 



(19) 



Note that T*,'^ = Tf^^ = ^ ^^^V and Tj^^^ = Tf" = 6.8 MeV to be consistent with eq.(14 ). e 
parametrizes the correction to the approximate identity of eq.(14). As shown in Appendix 
A , there is a model independent bound 

|e[ < 4% . (20) 

This correction introduces changes Li^a by up to ±2.0 fm^ but changes the other quantities 
by negligible amounts. 

Now it is clear that eqs. (17-19) are solvable. The three equations determine three quan- 
tities: Li^A and J dEpP„^/i,^, which are the magnitude of the total active neutrino 
flux, the axial two-body current, and the measured I'e/^x ratio, respectively. If there is no 
neutrino oscillation, then J dEpP^^/^^ — 1. 

For the experimental inputs, SNO has measured the ES rates, but the SK determination 
is more precise while in agreement with SNO. Therefore we use the SK measurements not 
only to provide a value for the integral above 6.8 MeV as shown above, but also to fix (and 
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TABLE I. Numbers of events reported by SNO for 306.4 live days for a kinetic energy threshold 
of 5 MeV. The uncertainties (except for the backgrounds) are statistical. 



Reaction 

Candidate Events 

Backgrounds 

Total Neutrino Events 

ES (from SK) 

Net NC + CC 

NC (CC shape unconstrained) 
NC (CC shape constrained) 



Events 
2928 
123 
2805 
258.3 
2546.7 
727 
576.5 



Uncertainty 

54T 
+21.6 -17.0 
+58.3 -56.7 

8.0 

59 

190 

49.5 



remove) the ES contribution to the total SNO rates above 5 MeV. With 1496 days of data, 
SK reports [27] the equivalent electron neutrino fluxes (= REs/'^e) for analysis thresholds 
of 6.5 and 5.0 MeV, 

$sk(6.5) = (2.3621^68) x lO' cm-^s"^ , (21) 
$sk(5.0) = (2.348^^66) X lO' cm-^s-^ , (22) 

respectively, where the statistical and systematic errors have been added in quadrature. Here 
and in the subsequent analysis in this paper, where asymmetric errors occur, we simply use 
the larger. We take the uncertainties in the two SK fluxes to be fully correlated. In view of 
the lack of significant threshold-energy dependence in the SK flux, we assume 

$sk(6.8) = $sk(6.5) . (23) 

SNO provides a model-independent value for the NC rate, obtained by separating the 
CC and NC parts of the signal by their different radial and sun-angle dependences in the 
detector. This result is independent of the energy spectrum of the CC events. Converting the 
flux reported into an equivalent number of events gives 727 ± 190 detected in the 306.4-day 
running period. 

In Table I the event rates needed for the model-independent analysis are summarized. 
The "true" number of ES events in the SNO data set is derived from $sk(5.0) and the SNO 
effective elastic scattering cross section with a 5.0 MeV threshold aelgg, 

CTelso = 1.10 X 10"^ cm^sT"^ , (24) 

derived from Ref. [1]. (The number is in excellent agreement with the 263.6 ± 26.4(stat.) 
obtained during the SNO signal extraction.) One could derive a value for the CC rate 
directly from the fifth and sixth lines of this table, but the two would be highly correlated. 
It is preferable to make use of expressions for NC + CC and NC because the summed rate 
is essentially free of correlation with the NC rate. So wc use 

Rnc + Rcc = Rtot - <I>sK(5.0)ae|5.o= (2546 ± 59) T"^ , 

Rj^c = (727 ± 190) T-i , (25) 
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where T=306.4 days. 

The averaged effective cross sections of SNO can be extracted from Ref. [1] 

a%c = 1-13 X 10-^ cm^sT-^ , 

a'^cc = 1-12 X 10"^ cm^sT-^ . (26) 

The effective cross sections are subject to uncertainty from a variety of sources, tabulated 
by SNO [1]. These include principally the energy scale, vertex-reconstruction accuracy, 
and (for the NC reaction) energy resolution and neutron capture efficiency. The sources of 
uncertainty produce in some cases correlated variations in the effective cross sections, which 
are explicitly accounted for in the analysis. 

In general we should have added 3% systematics for the NNLO EFT calculations of 
and o^c^ because for an EFT with a small expansion parameter Q ~ 1/3, 3% error 
is reasonable for a third order (NNLO) calculation. In the analysis of [2], however, a faster 
convergence is seen in four channels of (anti)ncutrino-deuteron scattering, such that 1-2% 
higher corrections also seems reasonable. Furthermore, it is conceivable that the higher order 
corrections can be absorbed in Li^a in low energy processes. One indication that this might 
happen is in the comparison with the potential model calculations. The potential model 
results have quite different systematics to those of EFT. The fact that NNLO EFT can fit 
four channels of (anti)neutrino-deuteron reaction results of [21] to within 1% [2] suggests 
that higher order effects can be absorbed in Li^a- Further investigation is still required to 
see whether the higher order effects shift Li^a approximately the same amount. For matrix 
elements with similar kinematics, this is likely to be true. In our case, we have CC and 
NC in approximately the same energy region. Thus we expect the higher order effects just 
shift Li A by a certain amount (~ +2 to +3 fm^, with the sign fixed by the fifth order 
proton-proton fusion calculation [16,3] as will be explained in more detail later) without 
introducing additional error to the and J dEpP^^i^^ determinations. 

For ES reactions with a 6.8-MeV threshold, the ratio of the neutral current and electron 
neutrino scattering cross sections is. 



0.153 , (27) 



with radiative corrections included. 

Now we have all the inputs required to solve eqs. (17-19). The full set of equations 
is nonlinear in Li^Ai but a linearized solution may be obtained by making a first-order 
expansion for gcc/9NC = 1 + (ctcc ~ ctArc)(Li^^/fm^ — 4.0) + . . . The term quadratic in 
(Li^^/fm^ — 4.0) is ~ 10~^ and can be neglected. Using this approximation, the solutions of 
eqs.(17-18) are: 



a 



NC 



X 



\[Rtot - *sK(5.0)ae[5o - Rnc] U ' ^\ + e) - $sK(6.8)a°,c + Rnc^^\ , (28) 



8 



TABLE II. Determinations of the NNLO (at renormalization scale mjr) from different 

processes. The higher order theoretical systcmatics are expected to be absorbed by shifting Li^a 
by ^ +2 to +3 fm^ thus is not included in this table. Note that the CC, NC &: ES combined 
analysis assumes the standard shape for the active neutrino flux. The tritium (3 decay analysis 
assumes the three-body current is negligible. The helioseismology analysis does not include the 
uncertainties from the solar model. The last two entries are theoretical determinations. EFT 
dimensional analysis gives \L\^a\ ^ 6 (fm^) which is denoted as [—6,6] as its expected range. 



Processes Li^a (fm^) References 

CC, NC & ES 4.0 ±6.3 [this work] 

Reactor V-d 3.6 ± 4.6 [3] 

Tritium (5 decay 4.2 ± 0.1 [29] (see also [16,3,30]) 

Helioseismology 4.8 ± 5.9 [31] 

Dimensional analysis ~ [—6, 6] [2] 

Potential model 4.0 [24] 



9nc(^ NC 



Inserting the experimental and theoretical quantities, 



/ 



Li,A = 4.0 ± 4.7(stat.) ± 4.5(syst.) fm^ , 

= (6.4 ± 1.4 ± 0.6) X 10^ cm-^s-^ , (31) 

dEpP,^,,^ = 0.25t°:J? ± 0.03 . 



The statistical errors in Li^a are dominated by RnCi and the systematic errors by $sk and 
by vertex reconstruction accuracy in SNO. 

A few comments can be made about the result we obtain in eq.(31). First, the only 
assumption we have made is that the active neutrino flux has the standard shape. The 
rest of the treatment is model independent in the sense that we have not assumed the size of 
Li A or assumed any neutrino oscillation scenarios. Second, our result on J dEpP^^/^^ can be 
used to constrain neutrino oscillation parameters. Third, the size of the active neutrino flux 
is consistent with the z/g flux of the standard solar model ^ssm = (S.OS^og}) x 10^ cm~^s~^. 
This sets a constraint on the oscillations between the active and sterile neutrinos. Fourth, 
the range of Li a we have obtained is consistent with the estimated value ~ 6 fm^ 

(at JI — ttItt) from dimensional analysis [2]. It is also consistent with the constraints from 
reactor-antineutrino deuteron breakup processes [3] , tritium beta decay [29] , helioseismology 
[31], and the latest improved potential model results [24] (corresponding to Li^a = 4.0 fm'^). 
The comparison of their corresponding NNLO Li^^'s is listed in Table 11. Here we have 
assumed that most of the higher order effects in EFT can be absorbed by Li^a, and the 
higher order theoretical systematics are therefore not included in the assigned error bars. 
We expect a +2 to +3 fm^ contribution to the effective value of Li^a from higher orders. The 
sign is fixed by an explicit fifth-order calculation of the proton-proton fusion at threshold 



9 



[16] which shows that Li^a shifts by +2 to +3 fm^ from the third order (NNLO) to the fifth 
order. Even though the tritium beta decay analysis assumes that the three-body current is 
neghgible and the hehoseismology analysis does not include the uncertainties from the solar 
model, it is still very encouraging that all the constraints agree with each other very well, 
given how different the physical systems are. 

It is likely in the future the error bar of Rnc could be reduced by a factor of 2. In that 
case, the error on Li^a would be reduced to 5 fm^. 

It is also interesting to reinvestigate the null hypothesis (specifically that all observed 
fluxes can be described consistently within the Standard Model of Particles and Fields) when 
Li^A is allowed to float. J dEpP^^/i,^ = 1 in the Standard Model, and thus the set of three 
equations (17-18) contain only two parameters. One finds that the set is inconsistent at 
4.3 (7. Alternatively, if one uses the experimental determination of Li^a from reactor data 
(Table II), the null hypothesis fails at 5.1 a (SNO only) or 5.3 a (SNO and SK). Thus, 
even if SNO were to place no rehance at all on the theoretical calculations of short-distance 
physics [21,24], it would still be true that the no-fiavor-conversion hypothesis is ruled out 
with high confidence. 

One might suspect that if the value of Li^a is taken from some other constraints, perhaps 
the shape assumption for the active neutrino flux can be removed. This question can be 
easily answered by inspecting the new set of equations 

Rnc — ^bo%c9nc j dEpNcPve-^v^ , 

Rcc = ^BO^cc9cc j dEpP,^^,^ , (32) 
Res = ^B^e ^ J dEp^^.P,^^,^ (1 + / ^^P^-e-.e 

where are the un-oscillated v^. flux and Pu^^ui is the probability distribution between 
the i/e — > i^i transition. If pjvc and p^^r satisfy the relation 

PArc|T*,^,=5MeV — P/i,T|T*^^=6.8MeV ) 

then one can determine j dEpNcPue^u^j J dEpP^^^^^ and $b provided Li^a is given. Un- 
fortunately, the above relation, which implies PAfc|r*f''^=5MeV — Pcc|T^'^=5MeV' <io^s not hold, 
as shown in Fig. 1(a) in Appendix A. 



IV. CONCLUSIONS 

We have analyzed the SNO and SK data on CC, NC and ES reactions to constrain 
the leading axial two-body current Li a- This two-body current contributes the biggest 
uncertainty in every low energy weak interaction deuteron breakup process, including SNO's 
CC and NC reactions. The only assumption made in this analysis is that the total flux of 
active neutrinos has the standard ^B spectral shape (but distortions in the electron neutrino 
spectrum are allowed). We have confirmed that SNO's conclusions about the inconsistency 
of the no-flavor-conversion hypothesis and the magnitude of the active solar neutrino flux 
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do not have significant theoretical model dependence. Our method has shown that SNO 
can be self-calibrated or be calibrated by SK with respect to theoretical uncertainties, and 
that the resulting calibration produces results in close accord with theoretical expectations. 
Alternatively, the purely experimental determination of Li a from reactor antineutrino data 
can be used to remove the dependence on theory, and SNO's conclusions are unaffected. 
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APPENDIX A: COMPUTATION OF THE NORMALIZED RESPONSE 
FUNCTIONS TO THE SPECTRUM 

In this Appendix, we define of the effective cross sections and then show the numerical 
results that support eqs. (12-14). 

For CC and ES reaction, the effective cross section a is related to the true cross section 
(7 through the relation 

^cciES) = Vcc(ES) f^^ dT j dT'rcciES)iT, T')^^^^ , (Al) 

'^CC(ES) 

where 77 is an experimental efficiency, T' is the true kinetic energy of the final state lepton 
and T is the electron energy recorded through the Cherenkov radiation of the electron, and 
T^^ is the detection threshold. If the resolution of the detector were perfect, the resolution 
function r{T,T') would be a delta function. For SNO, r{T,T') is a Gaussion function [1] 



rcc{ES){T,T') = ^ exp 

y Z-Kl\cC{ES) 



'/\2 



{T -T 



(A2) 



with resolution 



Acc = I -0.0684 + 0.331^/ ^ + 0.0425^ | MeV . (A3) 

For ES reactions, we have used the SK result rather than the SNO result for better statistics. 
The resolution for SK is [23] 



For NC reaction, the final state neutrons are thermalized then captured by deuterons to 
form tritons and photons in SNO's first phase running. The photons subsequently excite 
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electrons which produce Cherenkov radiation. Thus the SNO's NC events can be recorded 
as electron detections as well. However, the kinematic information of the final state neutrons 
is lost in thermalization. Thus the resolution function is monoenergetic 



rNc{T) 



1 



27rA 



■ exp 



NC 



NC) 



(A5) 



where T^c = 5.08 MeV and Aatc = l-H MeV. In contrast to eq.(A2), the effective NC 
differential cross section versus E is not distorted 



ctnc = Vncc^nc / dTrNc{T) 



Tth 



(A6) 



We now turn to the computation of the normalized response functions of the spectrum 
define in eq.(8). In Fig. 1(a), Pcc\T*h^=bMeV (sohd curve) and Pivc|T;»'p=5MeV (dashed curve) 
are shown as functions of E. The two curves are quite different, p^c is independent of Tj^^, 
according to eqs.(8) and (A6). In contrast, the peak of pcc can be shifted towards the high 
energy end by increasing Tqq. When T^'^ = 5MeV, pcc is close to a Gaussian function. 
Likewise, when T*'* = 6.8MeV, pe and p^^j. are adjusted to be close to Gaussians as well. 
Thus pcc-i Pe and p^,,- can be related. To see how different they are, it is convenient to define 
the following functions. 



^PCC = PCc|li,a=5 fmS — Pcc|li,a=0 ) 
5pe = Pe|T*''=6.8MeV ~ Pm.t lr*''=6.8MeV 
5pcC,ES = PCc\T^h 5MeV " Pe |T*''=6.8MeV 



(A7) 



5pcC) ^Pe and 5pcc.ES a-i'c shown as functions of E in Fig. 1 (b)-(d), respectively. 
To study the contributions of non-zero 5p, we will first prove an equality. Defining 



5p 



5p±\5p\ 



/± = y dE5p±P^^/^^ , 



then 



j dE5pP,^/,^^I+ + I^ . 

Since 0< P,,/,, {E)<1, 

< 7+ < J dE5p+ , 
J dE5p- <I-< . 
Because J dESp = 0, J dE5p± = ±| J dE \Sp\. Thus we find 



(A8) 



(A9) 



(AlO) 



/ 



dE5pPujy^ 



<- I dE\6p\ 



(All) 
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6 8 10 12 14 16 
E (MeV) 
(a) 



0.0004 



5 Pec 0002 



0.0002 




6Pe 



0.0004 
0.0002 



- 0.0002 




0.002 



5 Pec, 



ES 



0.002 

- 0.004 

- 0.006 




E (MeV) E (MeV) 

(c) (d) 
FIG. 1. (a) Pcc\T*^h^,=5MeV (solid curve) and pNc\Tj^^=5MeV (dashed curve) are shown as func- 
tions of E. (h)-(d) 5pcc, Spe and Spcc,ES defined in eq.(A7) are shown as functions of E, respec- 
tively. 



This model independent relation gives 



dE5pccPv 



dE5pEsPu,/u 



dE5p. 



CC,EsPu,/iy^ 



< 0.0013 



< 0.0012 



< 0.0105 



(A12) 



in comparison with J dEpPy^jy^ ~ 0.25. The first two inequality show that the Li^a de- 
pendence in pec and the difference between pe and are negligible compared with the 
correction from Spcc,ES- The Spcc,ES effect is parametrized by e in eq.(19) with 



e < 47o . 



(A13) 
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